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Abstract. We discuss the temporal evolution of the pair plasma created in Gamma-Ray Burst 
sources. A particular attention is paid to the relaxation of the plasma into thermal equilibrium. 
We also discuss the connection between the dynamics of expansion and the spatial geometry of the 
plasma. The role of the baryonic loading parameter is emphasized. 
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THERMALIZATION OF THE PAIR PLASMA 

In the recent publication [1] we considered the kinetic evolution of the pair plasma which 
is assumed to be formed in Gamma-Ray Burst (GRB) sources. Two alternative view- 
points exist in literature. In most papers, e.g. [2], it was assumed that the plasma quickly 
reaches thermal equilibrium, since the optical depth of the photons is much larger than 
one, and that it then expands adiabatically reaching a large Lorentz 7 factor. Cavallo and 
Rees [3], instead, suggested a number of different scenarios, assuming in all of them that 
the plasma cools down due to direct bremsstrahlung process until its temperature reaches 
0.511 MeV and that the electron-positron pairs disappear before starting to expand. In 
order to clarify this point we numerically integrated the relativistic Boltzmann equations, 
taking into account all 2-body and 3-body interactions between electrons, positrons and 
photons, assuming that the plasma is homogeneous and isotropic. We focused on the 
following range for the average energy E per particle 

0.1MeV<£ < lOMeV, (1) 

which is of interest for GRBs studies [4]. The results of our calculations are represented 
in Fig. 1. There, we also show how concentrations of pairs and photons evolve when 
inverse 3-body interactions are switched off artificially. 



FIREBALL VERSUS FIRESHELL 



GRBs are very different from the fireball phenomenon occurring e.g. in an atomic bomb 
explosion. This last one is characterized by a high temperature inside, a blast wave prop- 
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FIGURE 1. The density n of pairs and photons (left plot, with the total density in bold), their spectra 
dp/de (center plot) and their temperatures 9 —kT / inieC^) along with chemical potentials (p — ^/ [nieCp') 
(right plot) are shown, e is the energy of particle in units of electron rest mass energy nieC^. Two different 
initial conditions were considered: when only pairs are present with negligible amount of photons, and 
the opposite case (upper and lower figures respectively). In both cases the pair plasma relaxes to thermal 
equilibrium configuration on a timescale f,/, < 10^'^ s for our parameter range given in Eq.(l), i.e. much 
before it starts to expand on the timescale t < 10^^ s. We also show by dashed lines (on the upper 
left panel) the evolution of pairs and photons concentrations when the inverse 3 -body interactions are 
neglected. In this case we get a result similar to the one by Cavallo and Rees [3]. Further details are given 
inRef. [1]. 



agating into the surrounding dense atmosphere and a characteristic self-similar motion 
described by the Sedov solution. The latter can be easily obtained by considerations 
of dimensionality. In a GRB explosion, instead, it is created a dense, hot and optically 
thick plasma, made by electron-positron pairs and photons, which starts to expand adi- 
abatically being governed by equations of relativistic hydrodynamics. The surrounding 
medium is so rarefied that the motion of the plasma can be well represented by a rela- 
tivistic expansion into vacuum. In the laboratory frame (where the black hole is at rest), 
the pair plasma appears to expand as a narrow shell which keeps constant its thickness 
due to the Lorentz contraction [2]. Therefore we refer to it as a "fireshell" (see e.g. 
Refs. [5, 6] and references therein). However, we must here recall that photons emit- 
ted at the same time from the different regions of the fireshell surface reach a distant 
observer at different times, due to the ultrarelativistic expansion speed 7 > 10^ and the 
consequent Doppler contraction of the photons' arrival time at the detector. Therefore, 
a distant observer would not see the "fireshell" itself but the corresponding "EquiTem- 
poral Surfaces" (EQTSes), i.e. the surfaces locus of points source of photons reaching a 
distant observer at the same time [7, 8, 9, 10]. During the acceleration phase, occurring 
because of the large radiative pressure inside, the pair plasma engulfs a certain amount 
of baryonic matter, parameterized by 



5 = ^ < lo-^ 

Etot 



(2) 




Fireshell radius (cm) Fireshell radius (cm) 



FIGURE 2. The fireshell Lorentz factor (left panel) and temperature (right panel) are plotted versus 
the radius for selected values of the baryonic loading parameter. Filled circles show the transparency 
point. Squares show the moment of departure from thermal distributions of electron-positron pairs. The 
star denotes the moment when the temperature of the fireshell equals 511 keV. Calculations have been 
performed with the following initial conditions: total energy £,0, ~ \A1 ■ 10^^ erg, number density of 
pairs n± — 3.60 • lO"'^ cm^^, initial radius Rq = 1.74 • 10** cm. 



where Mb is the total mass of baryons and Etot is the total energy of plasma [11]. The 
whole system is now made by electrons, positrons, photons and protons and continues 
to accelerate until it becomes transparent, when the electron-positron pairs disappear 
and all photons escape ([2, 1 1], see Fig. 2). For simplicity, we assume a collision with a 
single baryonic matter shell placed between two given radii [11]. 

From hydrodynamic equations it is possible to show that the thickness of the fireshell 
indeed remains constant under the condition 7^ 1 [12]. This condition is approached 
soon after the beginning of expansion since the scaling law 7 r is valid for the 
radiation-dominated pair plasma. The above mentioned conclusion was obtained in [2] 
by analyzing different geometries of the plasma. Notice that the broadening of the shell 
does not take place if the strict inequality in Eq.(2) is satisfied [11, 13]. 

As one can see from Fig. 2, the influence of baryons on the dynamics of expansion 
is negligible for B < 10^. Besides, due to the high expansion rate, the concentration of 
electron-positron pairs is frozen out well before the moment of transparency [11]. 

In Fig. 3 we show the difference between numerical [2, 11] and analytical results 
[14, 15] for the Lorentz gamma factor reached at the transparency point [13]. While the 
results by [2] and [14] are similar for B > 10'^, marked differences arise in the range 
10^ < B < 10^, and the asymptotic values of gamma for B ^ are also different. 
Moreover, in the numerical treatment [11] this asymptotic behavior is reached for larger 
values of B, in disagreement with analytical expectations [14]. Thus, the acceleration of 
the fireshell for small B is larger if one accounts for pairs dynamics, which is neglected 
by the analytic treatments [13]. 
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FIGURE 3. The fireshell Lorentz gamma factor at the transparency point is plotted versus the B 
parameter. The thick line denotes numerical results [13], the plain line corresponds to analytical estimate 
from Shemi and Piran model [14], the dotted line denotes the asymptotic value y— B The dashed line 
shows results of Nakar, Piran and Sari [15]. 



CONCLUSIONS 

We conclude that the Cavallo and Rees fireball scenarios, a cavity filled of photons 
degraded to 0.5 IVIeV, never apply to GRBs. They systematically neglect the occur- 
rence of inverse 3-body interactions. If one takes into account both direct and inverse 
2-body and also 3-body interactions, one obtains thermal equilibrium on a much shorter 
timescale than the dynamical one. Thus, the initial state of the pair plasma is described 
uniquely by the common temperature of pairs and photons, while their chemical poten- 
tial is zero. Under these conditions the electron-positron-photon plasma is confined to a 
self-accelerating spherical shell, surrounding an almost empty cavity: the fireshell. 
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